Detection method of space domain maximum posteriori probability in a wireless communication system

ABSTRACT

A detection method of space domain maximum posteriori probability (MAP) in a wireless communication system is a kind of detection technique of space domain MAP in a multi-aerial wireless communication system, which makes the detection of iteration soft interference elimination for the interference signals in the other moments, and makes the MAP detection or the simplified MAP detection for the all signals in current time and space domain. It includes the steps of: taking a fading block as a unit; the received signals are matched and combined in time domain and space domain, and calculating an equivalent channel matrix; making the sequence detection, calculating estimation and variance of the signals before the ending of detecting internal iteration; when the internal iteration detection finishes, the result is outputted, or is outputted after hard decision.

CROSS-REFERENCE TO RELATED APPLICATION

This application is a National Phase Patent Application and claims the priority of International Application Number PCT/CN2007/000304, filed on Jan. 26, 2007, which claims priority of Chinese Patent Application Number 200610066285.4, filed on Mar. 31, 2006.

FIELD OF THE INVENTION

The present invention relates to a detection technique of space domain maximum a posteriori probability in a wireless communication system, and pertains to the technical field of high-speed wireless transmission.

BACKGROUND OF THE INVENTION

Employing multi-antenna transmitting and multi-antenna receiving technique in a wireless communication system can improve transmitting capacity of the communication system in many times theoretically. However, at the receiving end in a multi-antenna wireless communication system, signal interference in space domain (i.e., between antennas) exists. When signal transmission is carried out on a single carrier of wide band or multiple carriers of wide sub-bands, the wireless channel of each carrier becomes a frequency selective channel, i.e., inter-symbol interference in different times exists. Therefore, in a frequency selective channel environment, signal interference between different antennas, signal interference in different times, and Additive White Gaussian Noise (AWGN) exist at the receiving end of the multi-antenna system. In a multi-antenna system, the detector part of the receiver must recover the transmitted signals from the multi-antenna receiving signals superposed with signals from different transmitting antennae. The maximum a posteriori probability (MAP) detecting method for interference in time domain and space domain dimensions has unrealizable complexity in case there is a large number of transmitting antennae or multi-paths, and therefore can't be used in practical systems. In ordinary outdoor channel environments, usually direct paths exist, and correlation between antennae exists; for such channels, MMSE or ZF based detecting method can't achieve ideal performance.

In practical communication systems, usually error control encoding is utilized, in order to resist noise and interference; however, high-performance error control encoding methods usually employ soft decision decoding, i.e., the soft information of the bits must be provided, instead of hard decision decoding; therefore, the detector must provide soft information of the bits. At the receiving end, utilizing an iterative detection decoding receiver in which the detector works with the decoder in an iterative mode can greatly improve performance, when compared to a traditional receiver in which the detector and the decoder work with each other in a cascade mode. However, an iterative detection decoding receiver requires that the detector must take soft input and provide soft output, i.e., the detector must can not only soft output decision information to the decoder but also utilize the feedback result from the decoder as a priori information. It is an important task to seek for a soft input and soft output detector that has high performance but low complexity for multi-antenna wireless communication systems in a frequency selective channel environment with spatial correlation and direct path component, in order to support wide application of multi-antenna wireless communication systems.

Suppose the number of transmitting antennae is N, the number of receiving antennae is M, and the number of channel paths is L. In a frequency selective channel, suppose the complex baseband transmitting signal on antenna n in time k is s′_(n,k), the signal received on receiving antenna m is r′_(m,k), and set the signal response transmitted from transmitting antenna n to transmitting antenna m in path l is h_(m,n,l); then the relationship between transmitted baseband signals and received baseband signals is:

$\begin{matrix} {r_{m,k}^{\prime} = {{\sum\limits_{l = 1}^{L}{\sum\limits_{n = 1}^{N}{h_{{mn},l}s_{n,k}^{\prime}}}} + z_{m,k}}} & \lbrack 1\rbrack \end{matrix}$

It can be denoted in the following expression in matrix and vector form:

$\begin{matrix} {{r_{k}^{\prime} = {{\sum\limits_{l = 0}^{L - 1}{H_{l}s_{k - l}^{\prime}}} + z_{k}}}{{{wherein}:r_{k}^{\prime}} = \left\lbrack {r_{1,k}^{\prime},r_{2,k}^{\prime},\ldots\mspace{14mu},r_{M,k}^{\prime}} \right\rbrack^{T}},{s_{k}^{\prime} = \left\lbrack {s_{1,k}^{\prime},s_{2,k}^{\prime},\ldots\mspace{14mu},s_{N,k}^{\prime}} \right\rbrack^{T}},{z_{k} = \left\lbrack {z_{1,k},z_{2,k},\ldots\mspace{14mu},z_{M,k}} \right\rbrack^{T}},{H_{l} = {\begin{bmatrix} h_{11,l} & h_{12,l} & \ldots & h_{{1N},l} \\ h_{21,l} & h_{21,l} & \ldots & h_{{2N},l} \\ \vdots & \vdots & \ddots & \; \\ h_{{M\; 1},l} & \ldots & \; & h_{{MN},l} \end{bmatrix}.}}} & \lbrack 2\rbrack \end{matrix}$

SUMMARY OF THE INVENTION

Technical problem: the object of the present invention is to provide a space domain maximum a posteriori probability (MAP) detector in a multi-antenna wireless system in a frequency selective channel environment, which features with high performance and low complexity in implementation, and can meet the requirement of an iterative detection decoding receiver for soft input and soft output of detector.

Technical scheme: The detection method of space domain maximum a posteriori probability (MAP) in a wireless communication system in the present invention performs iterative soft interference cancellation detection for the interference signals in other moments, but performs maximum a posteriori probability (MAP) detection or simplified MAP detection for all signals in the current time in space domain, through the following steps:

-   1.1) Performing matched combination for received signals in time     domain and space domain on individual fading block basis, and     calculating an equivalent channel matrix; -   1.2) Performing sequential detection, carrying out interference     cancellation for the interference signals at other moments with     reference to the estimated value of the signal vector obtained in     step 1.1) to obtain received signals after interference     cancellation, and calculating an equivalent noise variance matrix,     and then performing MAP or simplified MAP detection with the     obtained signals after interference cancellation, equivalent channel     matrix, and equivalent noise variance matrix, to obtain the     detection result; -   1.3) Obtaining soft information of the bits from the detection     result obtained in step 1.2), Calculating estimated values and     variances of the signals before the internal iteration in the     detector ends, and returning to step 1.2) to perform next detection     cycle; outputting the result when the internal iteration in the     detector ends, or outputting the result after hard decision.     -   During the interference cancellation process in step 1.2), the         expectation of the interference signals is used as the estimated         value of the interference signals, the equivalent noise variance         matrix is calculated according to the channel matrix, variance         of noise, and variance of the interference signals; the         expectation and variance of the interference signals can be         calculated from the detection result in the last detection cycle         or from the soft information provided from the decoder.     -   Alternatively, the detecting method can be implemented in a         simplified manner through the following steps: -   3.1) Initialization: for each k, calculating:

${x_{k}^{\prime} = {\sum\limits_{l = 0}^{L - 1}{H_{l}^{H}r_{k + l}^{\prime}}}},$

-    and calculating

$\begin{matrix} {G^{\prime} = \begin{bmatrix} G_{- {({L - 1})}}^{\prime} & \ldots & G_{0}^{\prime} & \ldots & G_{L - 1}^{\prime} \end{bmatrix}} \\ {= {\begin{bmatrix} H_{0}^{H} & H_{1}^{H} & \ldots & H_{L - 1}^{H} \end{bmatrix}\begin{bmatrix} H_{L - 1} & \ldots & H_{0} & \; & \; \\ \; & \ddots & \vdots & \; & \; \\ \; & \; & \; & \; & \; \\ \; & \; & H_{L - 1} & \; & H_{0} \end{bmatrix}}} \end{matrix}$

-   G′₁ is a matrix obtained by removing G′₀ from G′, expanding the real     part and virtual Part of G′₀ and G′₁ to form new real matrixes,     i.e.,

${G_{0} = \begin{bmatrix} {{Re}\left( G_{0}^{\prime} \right)} & {- {{Im}\left( G_{0}^{\prime} \right)}} \\ {{Im}\left( G_{0}^{\prime} \right)} & {{Re}\left( G_{0}^{\prime} \right)} \end{bmatrix}},{{G_{l} = \begin{bmatrix} {{Re}\left( G_{l}^{\prime} \right)} & {- {{Im}\left( G_{l}^{\prime} \right)}} \\ {{Im}\left( G_{l}^{\prime} \right)} & {{Re}\left( G_{l}^{\prime} \right)} \end{bmatrix}};}$

-    expanding the real part and virtual part of x′_(k) as follows:

${x_{k} = \begin{bmatrix} {{Re}\left( x_{k}^{\prime} \right)} \\ {{Im}\left( x_{k}^{\prime} \right)} \end{bmatrix}};$

-    performing Cholesky factorization for G₀, to obtain G₀=U₁ ^(T)U₁,     wherein, U₁ is a upper triangular matrix composed of real numbers;     calculating {tilde over (G)}₁ and {tilde over (F)}₁, wherein, {tilde     over (G)}₁=U₁ ^(−T)G₁, each element of {tilde over (F)}₁ is the     square of the corresponding element in {tilde over (G)}₁; -   3.2) Sequential detection: performing the following steps in     sequence for k=1 . . . K: -   3.2.1) Interference cancellation: calculating {circumflex over     (x)}_(k)=x_(k)−G₁E[u_(1,k)], η_(k)={tilde over (F)}₁v_(k), wherein,     u_(1,k) is a vector composed of interference signals, v_(k) is a     vector composed of variances of interference signals; -   3.2.2) Calculating ŝ_(k)=G₀ ⁻¹{circumflex over (x)}_(k), for n=1, .     . . , 2N, calculating ρ_(n)=1/√{square root over (η_(k,n)+σ_(z) ²)},     and then multiplying the n^(th) row of U₁ by ρ_(n), to obtain a new     upper triangular matrix U; -   3.2.3) for all bits b_(i) corresponding to signal vector s_(k),     calculating:

${L_{D,E}\left( b_{i} \right)} = {{\min\limits_{\underset{{s_{k}:b_{i}} = {- 1}}{︸}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = 1}^{N}{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {\min\limits_{\underset{{s_{k}:b_{i}} = {+ 1}}{︸}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = 1}^{N}{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {L_{A}\left( b_{i} \right)}}$

-   3.2.4) If the current iteration is the last iteration cycle of     internal detection, outputting the result; otherwise reconstruct     estimated value and variance of the signals, and repeating step     3.2).

Functional effects: The detection method of space domain maximum a posteriori probability (MAP) in a multi-antenna wireless communication system in a frequency selective channel environment in the present invention has performance much superior to linear detection methods such as MMSE method when the channel condition is harsh, i.e., severe direct paths or strong correlation between antennae exist. But the complexity in implementation of the method is much lower than that of conventional MAP methods, especially in the cases that a large number of channel paths exist. At a lower complexity, the detection method settles the problem that MAP detectors can't be implemented in a multi-antenna system in a frequency selective channel environment due to extremely high complexity, and settles the problem that linear soft-input and soft-output detection methods have poor performance or event can't work normally under poor channel condition.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 are schematic diagrams of an encoding multi-antenna communication system and iterative receiver for space domain MAP detection.

FIG. 2 is a schematic diagram of soft-input and soft-output space domain MAP detector.

DETAILED DESCRIPTION OF THE EMBODIMENTS

The upper part of FIG. 1 illustrates a transmitter applicable to the communication system in the present invention, and the lower part of FIG. 1 illustrates an iterative receiver applicable to the communication system in the present invention. FIG. 2 is a schematic diagram of soft-input and soft-output space domain MAP detector. The technical scheme mainly comprises four steps: initialization, interference cancellation, MAP detection, and calculation of symbol probability and statistical values. Hereunder the four steps will be detailed, and the workflow of the detector will be discussed finally.

1. Initialization:

The detector works on individual block basis, i.e., in a block, the channel is deemed as time-independent. Suppose the block length is K. The initialization comprises the following steps in sequence:

-   (1.1) Time-space combination: for each k=1, . . . , K, calculate:

$x_{k}^{\prime} = {\sum\limits_{l = 0}^{L - 1}{H_{l}^{H}r_{k + l}^{\prime}}}$

-   (1.2) Calculate an equivalent channel coefficient matrix:

$\begin{matrix} {G^{\prime} = \left\lbrack {G_{- {({L - 1})}}^{\prime}\mspace{14mu}\ldots\mspace{14mu} G_{0}^{\prime}\mspace{14mu}\ldots\mspace{14mu} G_{L - 1}^{\prime}} \right\rbrack} \\ {= {\left\lbrack {H_{0}^{H}H_{1}^{H}\mspace{14mu}\ldots\mspace{14mu} H_{L - 1}^{H}} \right\rbrack\begin{bmatrix} H_{L - 1} & \ldots & H_{0} & \; \\ \; & \ddots & \vdots & \; \\ \; & \; & \; & \; \\ \; & \; & H_{L - 1} & H_{0} \end{bmatrix}}} \end{matrix}$

Let G′₁ is a matrix obtained by removing G′₀ from G′,

-   (1.3) Expansion of complex signals into real signals: expand the     real part and virtual part of G′₀, G′₁, and x′_(k), i.e.:

${G_{0} = \begin{bmatrix} \begin{matrix} {{Re}\left( G_{0}^{\prime} \right)} & {- {{Im}\left( G_{0}^{\prime} \right)}} \end{matrix} \\ \begin{matrix} {{Im}\left( G_{0}^{\prime} \right)} & {{Re}\left( G_{0}^{\prime} \right)} \end{matrix} \end{bmatrix}},{G_{I} = \begin{bmatrix} \begin{matrix} {{Re}\left( G_{I}^{\prime} \right)} & {- {{Im}\left( G_{I}^{\prime} \right)}} \end{matrix} \\ \begin{matrix} {{Im}\left( G_{I}^{\prime} \right)} & {{Re}\left( G_{I}^{\prime} \right)} \end{matrix} \end{bmatrix}},{x_{k} = \begin{bmatrix} {{Re}\left( x^{\prime} \right)} \\ {{Im}\left( x^{\prime} \right)} \end{bmatrix}},$

-   (1.4) Cholesky factorization: G₀=U₁ ^(T)U₁, wherein, U₁ is a upper     triangular matrix composed of real numbers, -   (1.5) Calculate {tilde over (G)}₁ and {tilde over (F)}₁, wherein,     {tilde over (G)}₁=U₁ ^(−T)G₁, each element of {tilde over (F)}₁ is     the square of the corresponding element in {tilde over (G)}₁;     2. Interference Cancellation:     -   In this step, interference cancellation is carried out for the         signals at other moments, and then the variance of interference         noises is calculated; next, a new upper triangular matrix is         calculated with the variance matrix, and the zero-forcing (ZF)         solution of the signals is calculated. The detail steps as         follows: -   (2.1) Interference cancellation and calculation of variance of     interference noise: calculate {circumflex over     (x)}_(k)=x_(k)−G₁E[u_(1,k)], wherein, u_(1,k)=[s_(k−L+1) ^(T), . . .     , s_(k−1) ^(T), s_(k+1) ^(T), . . . , s_(k+L−1) ^(T)]^(T). Calculate     η_(k)={tilde over (F)}₁v_(k), wherein, v_(k) is a column vector, and     each element in v_(k) is equal to the variance of the corresponding     element in u_(1,k). -   (2.2) Calculate a new upper triangular matrix: for n=1, . . . , 2N,     calculate ρ_(n)=1/√{square root over (η_(k,n)+σ_(z) ²)}, and then     multiply the n^(th) row of U₁ by ρ_(n), to obtain a new upper     triangular matrix U; -   (2.3) Calculate the ZF solution of the signals: ŝ_(k)=G₀     ⁻¹{circumflex over (x)}_(k).     3. MAP Detection:

In this step, the likelihood ratios (i.e., external detection information) of specific bits are calculated according to the upper triangular matrix U, ZF solution ŝ_(k) and a priori probability P(s_(n,k)).

$\begin{matrix} {{L_{D,E}\left( b_{i} \right)} = {{\min\limits_{\underset{{s_{k}\text{:}b_{i}} = {- 1}}{︸}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = 1}^{2\; N}{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {\min\limits_{\underset{{s_{k}\text{:}b_{i}} = {+ 1}}{︸}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = 1}^{2\; N}{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {L_{A}\left( b_{i} \right)}}} & (3) \end{matrix}$

If expression [3] is calculated with an exhaustion method, the complexity of computation will increase exponentially as N increases. In the present invention, a search method with lower complexity is used, as follows:

-   (3.1) Search for a vector b, in which each element is ±1, and the     real signal vector s obtained by mapping b will make the value of

${{U\left( {{\hat{s}}_{k} - s} \right)}}^{2} - {\sum\limits_{n = 1}^{2\; N}{\log\;{P\left( s_{n} \right)}}}$

-    minimum. This step can be implemented through the following     process: -   (3.1.1) For i=1, . . . , 2N, j=i+1, . . . , 2N, q_(ii)=u_(ii) ²,     q_(ij)=u_(ij)/u_(ii); let n=2N; C_(opt)=+∞; C_(cur)=0; θ_(n)=0;     t_(n)=1;

For each possible value α_(l) of S_(n), calculate μ_(n,l)=q_(nn)(α_(l)−ŝ_(n))²−log P(s_(n)=α_(l)), sort μ_(n,l) in ascending order, to obtain a group of sequence numbers d(l) that makes μ_(n,d(1))<μ_(n,d(2))< . . . μ_(n,d(M) _(C) ₎, wherein, M_(C) is the number of bits corresponding to each real sign.

-   (3.1.2) Calculate C_(cur)=θ_(n)+μ_(n,d(t) _(n) ₎, -   If C_(cur)<C_(opt), and n>1, then go to 3.1.3) -   If C_(cur)<C_(opt), and n=1, then go to 3.1.4) -   If C_(cur)≧C_(opt), and n=2N, then go to 3.1.6) -   If C_(cur)≧C_(opt), and n<2N, then go to 3.1.5) -   (3.1.3) n=n−1, θ_(n)=C_(cur); t_(n)=1; for each possible signal     value α_(l) of s_(n), calculate:

$\mu_{n,l} = {{q_{nn}\left( {\alpha_{l} - {\hat{s}}_{n} + {\sum\limits_{j = {n + 1}}^{2\; N}{q_{nj}\left( {s_{j} - {\hat{s}}_{j}} \right)}}} \right)}^{2} - {\log\;{P\left( {s_{n} = \alpha_{l}} \right)}}}$ Sort μ_(n,l) in ascending order to obtain a group of sequence numbers d(l), which makes μ_(n,d(1))<μ_(n,d(2))< . . . μ_(n,d(M) _(C) ₎; go to 3.1.2;

-   (3.1.4) s_(map)=s; C_(opt)=C_(cur); n=n+1; t_(n)=t_(n)+1; if     t_(n)>M_(C), μ_(n,d(t) _(n) ₎=+∞; go to 3.1.2; -   (3.1.5) n=n+1; t_(n)=t_(n)+1; if t_(n)>M_(C), μ_(n,d(t) _(n) ₎=+∞;     go to 3.1.2; -   (3.1.6) Map s_(map) reversely to b_(map); output the resulting     b_(map) and C_(map)=C_(opt). -   (3.2) For i=1 . . . 2N, search for a bit vector {circumflex over     (b)} respectively, wherein, the i^(th) bit vector {circumflex over     (b)}_(i) meets {circumflex over (b)}_(i)=−b_(map,i) and ensures the     real signal vector s obtained by mapping {circumflex over (b)} makes     the value of

$C_{l} = {{{U\left( {{\hat{s}}_{k} - s} \right)}}^{2} - {\sum\limits_{n = 1}^{2\; N}{\log\;{P\left( s_{n} \right)}}}}$

-   -    minimum; output the minimum C_(i). This step can be implemented         as follows:

-   (3.2.1) Let n=2N; C_(opt)=+∞; C_(cur)=0; θ_(n)=0; t_(n)=1;

For each possible value α_(l) of s_(n) that meets {circumflex over (b)}_(i)=−b_(map,i), calculate μ_(n,l)=q_(nn)(α_(l)−ŝ_(n))²−log P(s_(n)=α_(l)) and sort μ_(n,l) in ascending order, to obtain a group of sequence numbers d(l) that makes μ_(n,d(1))<μ_(n,d(2))< . . . μ_(n,d(M) _(C) _(′)); if s_(n) contains {circumflex over (b)}_(i), then M_(C)′=M_(C)/2; if s_(n) doesn't contain {circumflex over (b)}_(i), then M_(C)′=M_(C).

-   (3.2.2) Calculate C_(cur)=θ_(n)+μ_(n,d(t) _(n) ₎, -   If C_(cur)<C_(opt), and n=1, then go to 3.2.3) -   If C_(cur)<C_(opt), and n=1, then go to 3.2.4) -   If C_(cur)≦C_(opt), and n=2N, then go to 3.2.6) -   If C_(cur)≦C_(opt), and n=2N, then go to 3.2.5) -   (3.2.3) n=n−1; θ_(n)=C_(cur); t_(n)=1; for each possible signal     value α₁ of s_(n) that meets {circumflex over (b)}_(i)=−b_(map),     calculate

$\mu_{n,l} = {{q_{nn}\left( {\alpha_{l} - {\hat{s}}_{n} + {\sum\limits_{i = {n + 1}}^{2\; N}{q_{nj}\left( {s_{j} - {\hat{s}}_{j}} \right)}}} \right)}^{2} - {\log\;{P\left( {s_{n} = \alpha_{l}} \right)}}}$ Sort μ_(n,l) in ascending order, to obtain a group of sequence numbers d(l) that makes μ_(n,d(1))>μ_(n,d(2))> . . . μ_(n,d(M) _(C) ′); go to (3.2.2).

-   (3.2.4) s_(map)=s; C_(opt)=C_(cur); n=n+1; t_(n)=t_(n)+1; if     t_(n)>M_(c)′; μ_(n,d(t) _(n)) =+∞; go to (3.2.2). -   (3.2.5) n=n+1; t_(n)=t_(n)+1; if t_(n)>M_(c)′, μ_(n,d(t) _(n)) =+∞;     go to (3.2.2); -   (3.2.6) C_(i)=C_(opt). -   (3.3) For each bit b_(i), calculate L_(D,E)(b_(i))={circumflex over     (b)}_(map,l)(C_(l)−C_(map))−L_(A)(b_(l)).

In the MAP detection step, the likelihood ratio of bit can be obtained with above method. When the detection iteration for a block is completed, the likelihood ratio can be outputted to the cecoder for soft decision decoding; before the detection iteration for a block is completed, the likelihood ratio can be used to reconstruct average value and variance of the signals; then, the interference cancellation step and the MAP detection step can be repeated.

(4.) Calculation of symbol probability and statistical values:

In the internal iteration process in the detector or the iteration process between the detector and the decoder, the likelihood ratios of bits must be converted to statistical values of symbols, i.e., expectation and variance, so as to perform interference cancellation. In addition, in the MAP detection procedure, a signal equal to the logarithm of probability of the symbols in the symbol set is required; the probability can be reconstructed from the soft information of bits from the decoder.

(4.1) Calculate symbol probability from soft information of bits

Suppose a bit sequence b_(n,k,j), j=1, . . . , M_(C) is mapped to symbol α, then the symbol probability can be calculated from soft information with the following expression:

${P\left( {s_{n,k} = \alpha} \right)} = {\prod\limits_{j = 1}^{M_{c}}\;{P\left( b_{n,k,j} \right)}}$ (4.2) Calculate statistical values of symbols with symbol probability:

The expectation value can be calculated with the following expression:

${E\left\lbrack s_{n,k} \right\rbrack} = {\sum\limits_{\alpha}{\alpha\;{P\left( {s_{n,k} = \alpha} \right)}}}$ The variance value can be calculated with the following expression:

${{cov}\left\lbrack s_{n,k} \right\rbrack} = {{\sum\limits_{\alpha}{{\alpha }^{2}\;{P\left( {s_{n,k} = \alpha} \right)}}} - {E\left\lbrack s_{n,k} \right\rbrack}^{2}}$ 5. Workflow of detector

In the initial detection and decoding cycle of the iterative detection and encoding receiver, since no soft information from the decoder is available, the detector performs iteration internally first. In the subsequent detection and decoding process, the decoder provides soft information, the detector doesn't any longer perform internal iteration; instead, the iteration is performed between the detector and the decoder. The present invention provides a soft-input and soft-output detection method of space domain maximum a posteriori probability in an multi-antenna wireless communication system in a frequency selective channel environment, which features with lower complexity and robust performance. The implementation is as follows:

-   1) Determine the system parameters, such as the number of     transmitting antennae, number of receiving antennaie, maximum number     of multi-paths of channel, and error control encoding mode, etc. In     this embodiment, the number of transmitting antennaie is 4, the     number of receiving antennae is 4, the number of multi-paths of     channel is 6, and the error control encoding mode is Turbo coding     with ½ code rate. -   2) Determine the number of self-iterations of the detector according     to the requirement of the receiving end for complexity and     performance: if the detector is iterative, the number of iterations     between the detector and the decoder must be determined further. In     this embodiment, the number of self-iterations of the detector     during initial detection decoding is 2; the number of iterations     between the detector and the decoder is 3. -   3) Design a soft-inpout and soft-out detection method according to     technical scheme 1-5, and implement through the following steps: -   3.1) Take the received signals in a fadin gblock as the object, and     perform initialization with the method described above. -   3.2) Take the initialized signals in a fading block as the object,     and perform two times iterative sequential detection in the     detector; the sequential detection comprises three steps:     interference cancellation, space domain MAP detection, and     calculation of statistical values of the symbols, which are     implemented through the steps of interference cancellation space     domain MAP detection, and calculation of symbol probability and     statistical values of the symbols as described above. -   3.3) Perform three times iterative sequential detection between the     detector and the decoder, the decoder can employ any soft-input and     soft-output decoding method; the sequential detection comprises     three steps: interference cancellation, space domain MAP detection,     and calculation of statistical values of the symbols, which are     implemented through the steps of interference cancellation, space     domain MAP detection, and calculation of symbol probability and     statistical values of the symbols as described above. -   3.4) Output the final result. 

1. A detection method of space domain maximum a posteriori probability (MAP) in a wireless communication system comprising a detector configured to perform an internal iteration, which is characterized in that it performs iterative soft interference cancellation detection for interference signals in other moments, but performs maximum a posteriori probability (MAP) detection or simplified MAP detection for all signals in a current time in space domain, through the following steps: 1.1) performing matched combination of received signals in time domain and space domain on individual fading block basis to compute a signal vector, and calculating an equivalent channel matrix; 1.2) perfoiming sequential detection, carrying out interference cancellation for the interference signals at other moments with reference to an estimated value of the signal vector obtained in step 1.1) to obtain received signals after interference cancellation, and calculating an equivalent noise variance matrix, and then performing the MAP or simplified MAP detection with the obtained received signals after interference cancellation, the equivalent channel matrix, and the equivalent noise variance matrix, to obtain a detection result comprising soft information of bits; and 1.3) obtaining the soft information of the bits from the detection result obtained in step 1.2), calculating an estimated value and a variance of the signals before the end of the internal iteration in the detector, and returning to step 1.2) to perform a next detection cycle; outputting a result at the end of the internal iteration in the detector, or outputting the result after a hard decision.
 2. The detection method of space domain maximum a posteriori probability (MAP) in a wireless communication system according to claim 1, wherein, during the interference cancellation process in step 1.2), the expectation of the interference signals is used as an estimated value of the interference signals, the equivalent noise variance matrix is calculated according to the equivalent channel matrix, the variance of noise, and the variance of the interference signals; and wherein the expectation and variance of the interference signals are calculated from the detection result in the last detection cycle or from the soft information.
 3. The detection method of space domain maximum a posteriori probability (MAP) in a wireless communication system according to claim 1, 3.1) wherein the calculating the equivalent channel matrix G′ comprises, k=1 . . . K, calculating: ${{x^{\prime}}_{k} = {\sum\limits_{l = 0}^{L - 1}\;{H_{l}^{H}r_{k + l}^{\prime}}}},$  and calculating: $\begin{matrix} {G^{\prime} = \begin{bmatrix} G_{- {({L - 1})}}^{\prime} & \ldots & G_{0}^{\prime} & \ldots & G_{L - 1}^{\prime} \end{bmatrix}} \\ {= {\begin{bmatrix} H_{0}^{H} & H_{1}^{H} & \ldots & H_{L - 1}^{H} \end{bmatrix}\begin{bmatrix} H_{L - 1} & \ldots & H_{0} & \; & \; \\ \; & \ddots & \vdots & \; & \; \\ \; & \; & H_{L - 1} & \; & H_{0} \end{bmatrix}}} \end{matrix}$  wherein G′₁ is a matrix obtained by removing G′₀ from G′, expanding the real part and imaginary part of G′₀ and G′₁ to form new real matrixes, that is, ${G_{0} = \begin{bmatrix} {{Re}\left( G_{0}^{\prime} \right)} & {- {{Im}\left( G_{0}^{\prime} \right)}} \\ {{Im}\left( G_{0}^{\prime} \right)} & {{Re}\left( G_{0}^{\prime} \right)} \end{bmatrix}},{{G_{l} = \begin{bmatrix} {{Re}\left( G_{l}^{\prime} \right)} & {- {{Im}\left( G_{l}^{\prime} \right)}} \\ {{Im}\left( G_{l}^{\prime} \right)} & {{Re}\left( G_{l}^{\prime} \right)} \end{bmatrix}};}$  expanding the real part and imaginary part of x′_(k) as follows: $x_{k} = \left\lbrack \begin{matrix} {{Re}\left( x_{k}^{\prime} \right)} \\ {{{Im}\left( x_{k}^{\prime} \right)};} \end{matrix} \right.$  performing Cholesky factorization for G₀, to obtain G₀=U₁ ^(T)U₁, wherein, U₁ is an upper triangular matrix composed of real numbers; calculating {tilde over (G)}₁ and {tilde over (F)}₁, wherein, {tilde over (G)}₁=U₁ ^(−T)G₁, each element of {tilde over (F)}₁ is the square of the corresponding element in {tilde over (G)}₁; and 3.2) wherein the performing the sequential detection is performed for k=1 . . . K: 3.2.1) wherein the interference cancellation comprises calculating {circumflex over (x)}_(k)=x_(k)−G₁E[u_(i,k], η) _(k)={tilde under (F)}₁v_(k), wherein, u_(l,k) is a vector composed of the interference signals, v_(k) is a vector composed of variances of the interference signals; 3.2.2) wherein the calculating the equivalent noise matrix comprises calculating ŝ_(k)=G₀ ⁻¹{circumflex over (x)}_(k), for n=1, . . . , 2N, calculating ρ_(n)=1/√{square root over (η_(k,n)+σ_(z) ²)}, and then multiplying the n^(th) row of U₁ by ρ_(n), to obtain a new upper triangular matrix U; 3.2.3) wherein the performing the MAP or simplified MAP detection comprises, for all bits b_(i), corresponding to the signal vector s_(k), calculating: ${{L_{D,E}\left( b_{i} \right)} = {{\frac{\min}{{s_{k}:b_{i}} = {- 1}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = l}^{2N}\;{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {\frac{\min}{{s_{k}:b_{i}} = {+ 1}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = l}^{2N}\;{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {L_{A}\left( b_{i} \right)}}};$  and 3.2.4) if the current iteration is the last iteration cycle of internal detection, outputting the result; otherwise reconstructing estimated value and variance of signals, and repeating the sequential detection, wherein K is the block length, L is the number of channel paths, N is the number of transmitting antennae, M is the number of receiving antennae, s′_(n,k) is the complex baseband transmitting signal on antenna n in a k-th block, r′_(k) is the signal received on a receiving antenna in a k-th block, σ_(z) ² is the variance of the noise in the signal received r′_(k); h_(mn,l) is the signal response transmitted from a transmitting antenna n to a receiving antenna m over path I, H_(l) is the channel matrix for the l-th tap $\begin{bmatrix} H_{11,l} & h_{12,l} & \cdots & h_{{1N},l} \\ H_{21,l} & h_{21,l} & \cdots & H_{{2N},l} \\ \vdots & \vdots & ⋰ & \; \\ h_{{M1},l} & \cdots & \; & H_{{M\; N},l} \end{bmatrix},$  and H_(l) ^(H) is the Hermitian transpose of H_(l).
 4. The detection method of space domain maximum a posteriori probability (MAP) in a wireless communication system according to claim 2: 4.1) wherein the calculating the equivalent channel matrix G′ comprises, for each k, calculating: ${{x^{\prime}}_{k} = {\sum\limits_{l = 0}^{L - 1}\;{H_{l}^{H}r_{k + l}^{\prime}}}},$  and calculating: $\begin{matrix} {G^{\prime} = \begin{bmatrix} G_{- {({L - 1})}}^{\prime} & \ldots & G_{0}^{\prime} & \ldots & G_{L - 1}^{\prime} \end{bmatrix}} \\ {= {\begin{bmatrix} H_{0}^{H} & H_{1}^{H} & \ldots & H_{L - 1}^{H} \end{bmatrix}\begin{bmatrix} H_{L - 1} & \ldots & H_{0} & \; & \; \\ \; & \ddots & \vdots & \; & \; \\ \; & \; & H_{L - 1} & \; & H_{0} \end{bmatrix}}} \end{matrix}$  wherein G′₁ is a matrix obtained by removing G′₀ from G′; expanding the real part and imaginary part of G′₀ and G′₁ to form new real matrixes, that is, ${G_{0} = \begin{bmatrix} {{Re}\left( G_{0}^{\prime} \right)} & {- {{Im}\left( G_{0}^{\prime} \right)}} \\ {{Im}\left( G_{0}^{\prime} \right)} & {{Re}\left( G_{0}^{\prime} \right)} \end{bmatrix}},{{G_{l} = \begin{bmatrix} {{Re}\left( G_{l}^{\prime} \right)} & {- {{Im}\left( G_{l}^{\prime} \right)}} \\ {{Im}\left( G_{l}^{\prime} \right)} & {{Re}\left( G_{l}^{\prime} \right)} \end{bmatrix}};}$  expanding the real part and imaginary part of x′_(k) as follows: $x_{k} = \left\lbrack \begin{matrix} {{Re}\left( x_{k}^{\prime} \right)} \\ {{{Im}\left( x_{k}^{\prime} \right)};} \end{matrix} \right.$  performing Cholesky factorization for G₀, to obtain G₀=U₁ ^(T)U₁, wherein, U₁ is an upper triangular matrix composed of real numbers; calculating {tilde over (G)}₁ and {tilde over (F)}₁, wherein, {tilde over (G)}₁=U₁ ^(−T)G₁, each element of {tilde over (F)}₁ is the square of the corresponding element in {tilde over (G)}₁; and 4.2) wherein the performing the sequential detection is performed for k=1 . . . K: 4.2.1) wherein the interference cancellation comprises calculating {circumflex over (x)}_(k)=x_(k)−G₁E[u_(i,k)], η_(k)={tilde over (F)}₁v_(k), wherein, u_(l,k) is a vector composed of the interference signals, v_(k) is a vector composed of variances of the interference signals; 4.2.2) wherein the calculating the equivalent noise matrix comprises calculating ŝ_(k)=G₀ ⁻¹{circumflex over (x)}_(k), for n=1, . . . , 2N, calculating ρ_(n)=1/√{square root over (η_(k,n)+σ_(z) ²)}, and then multiplying the n^(th) row of U₁ by ρ_(n), to obtain a new upper triangular matrix U; 4.2.3) wherein the performing the MAP or simplified MAP detection comprises, for all bits b_(i) corresponding to the signal vector s_(k), calculating: ${{L_{D,E}\left( b_{i} \right)} = {{\frac{\min}{{s_{k}:b_{i}} = {- 1}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = l}^{2N}\;{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {\frac{\min}{{s_{k}:b_{i}} = {+ 1}}\left\{ {{{U\left( {{\hat{s}}_{k} - s_{k}} \right)}}^{2} - {\sum\limits_{n = l}^{2N}\;{\log\;{P\left( s_{n,k} \right)}}}} \right\}} - {L_{A}\left( b_{i} \right)}}};$  and 4.2.4) If the current iteration is the last iteration cycle of internal detection, outputting the result; otherwise reconstructing estimated value and variance of signals, and repeating the sequential detection, wherein K is the block length, L is the number of channel paths, N is the number of transmitting antennae, M is the number of receiving antennae, s′_(n,k) is the complex baseband transmitting signal on antenna n in a k-th block, r′_(k) is the signal received on a receiving antenna in a k-th block, σ_(z) ² is the variance of the noise in the signal received r′_(k), h_(mn,l) is the signal response transmitted from a transmitting antenna n to a receiving channel matrix for the l-th tap $\begin{bmatrix} H_{11,l} & h_{12,l} & \cdots & h_{{1N},l} \\ H_{21,l} & h_{21,l} & \cdots & H_{{2N},l} \\ \vdots & \vdots & ⋰ & \; \\ h_{{M1},l} & \cdots & \; & H_{{M\; N},l} \end{bmatrix},$  and H_(l) ^(H) is the Hermitian transpose of H_(l). 